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LAUNCHER  DYNAMICS  STUDY 
Part  R 


Approximate  Formulas  for  Pitch  Velocity  at  Tip-off 


ABSTRACT: 

The  preliminary  design  or  feasibility  study  of  a  straight  rail  launching 
system  requires,  in  addition  to  layouts  and  material  composition,  estimates  of 
the  accuracy  of  the  system.  Possible  sources  of  inaccuracy,  which  is  defined  to 
be  nonreproducable  tip-off  velocity,  are  thrust  misalignment,  rail  curvature 
arising  from  manufacturing  processes  or  uneven  heating  due  to  the  sun,  shaking 
forces  of  unknown  phase  which  might  exist  in  SOSR  designs,  or  excessive  amounts 
of  friction  between  shoes  and  rail. 

Our  purpose  in  this  report  is  to  present  techniques  of  analysis  for  a 
straight  rail  system  with  a  view  toward  obtaining  estimates  of  changes  in  pitch 
velocity  of  the  rocket  at  tip-off  as  might  exist  because  of  the  presence  of  these 
various  disturbing  agencies. 
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Additional  symbols  are  defined  in  the  body  of  the  text. 


I.  INTRODUCTION 


The  preliminary  design  or  feasibility  study  of  a  straight  rail  launching 
system  requires,  in  addition  to  layouts  and  material  composition,  estimates  of 
its  accuracy.  Possible  sources  of  inaccuracy,  which  is  defined  to  be  nonreproducabTe 
tip-off  pitch  velocity,  are  thrust  misalignment,  rail  curvature  arising  from 
manufacturing  processes  or  uneven  heating  due  to  the  sun,  shaking  forces  which 
mig^t  exist  in  SOSR  designs,  or  excessive  amounts  of  friction  between  shoes  and 
rail. 

Our  purpose  in  this  report  is  to  present  techniques  of  analysis  for  a 
straight  rail  system  with  a  view  toward  obtaining  estimates  of  changes  in  pitch 
velocity  of  the  missile  at  tip-off  as  mi^t  exist  because  of  the  presence  of  these 
various  disturbing  agencies.  Since  these  results  are  to  be  used  for  feasibility 
studies  or  the  like,  methods  will  be  developed  so  that  these  quantities  can  be 
obtained  by  means  of  slide  rule  or  desk  machine  calculations  rather  than  by 
means  of  analogue  or  digital  computing  devices.  For  the  detailed  analyses 
attendent  to  final  design  proposals,  numerical  investigations  of  the  type  carried 
out  in  many  of  the  previous  technical  reports  of  this  project  are  necessary. 

All  of  these  methods  are  associated  with  a  launcher  whose  gross  con¬ 
figuration  can  be  accurately  represented  by  the  two  degree  of  freedom  model 
which  appeared  in  TAM  Report  No.  552,  University  of  Illinois,  Urbana,  Illinois, 

1960.  As  shown  in  Fig.  1  ,  this  model  consists  of  a  rigid  L-frame  supported 
at  a  single  point  O  .  Since  this  is  a  pivot  point,  horizontal  and  vertical  motion 
is  prevented  and  so  only  rotation  is  possible.  Stability  is  provided  by  some 
equivalent  helical  spring  which  represents  the  overall  flexibility  of  the  structure. 

The  launching  rail  AB  may  have  some  small  curvature  as  indicated. 


FIGURE  1 


The  missile  is  represented  by  two  masses  rigidly  attached  to  two  sleeves 
that  slide  along  the  arm  of  the  launching  frame.  These  masses  are  the  source 
of  rotational  inertia  for  the  missile.  (We  only  consider  systems  in  which  the 
shoes  of  the  rocket  leave  the  rail  simultaneously.)  Finally  the  thrust  force  is 
denoted  by  T  (t). 

This  dynamical  system  has  two  degrees  of  freedom.  The  coordinates 
used  to  describe  its  configuration  at  time  t  are  the  angular  rotation  0  (t)  of 
its  arm  as  measured  from  equilibrium  and  (t)  ,  the  distance  the  mass  center 
of  the  missile  has  moved  up  the  rail.  As  in  all  previous  work,  it  is  again  assumed 
that  the  angle  6  is  small  so  that  the  usual  approximations  based  on  linearization 
are  permissible. 


3. 


The  differential  equations  which  describe  the  motion  of  this  system  are 
very  complex.  We  may;  however,  in  certain  instances,  bypass  complete  integra¬ 
tion  of  the  equations  of  motion  of  any  system.  Our  efforts  here  are  based  on  the 
reduction  of  the  differential  equations  of  motion  to  integral  equations  and  a  sub¬ 
sequent  solution  by  iteration.  The  necessary  number  of  steps  in  the  Iteration 
depends  on  the  numerical  values  of  a  parameter  which  is  the  ratio  launch  time 
versus  natural  period  of  the  system.  In  some  instances  we  can  employ  the 
principle  of  Impulse  and  momentum  to  interpret  these  approximations. 

Equations  of  Motion  for  the  Model 

We  begin  with  a  discussion  of  the  kinematics  of  the  model  shown  schemati¬ 
cally  in  Fig.  2  .  The  missile  Is  represented  by  the  bent  ABCD  and  the  launcher 
by  the  rigid  frame  ORS  .  Three  coordinate  systems  are  introduced.  The  first 
of  these  labelled  XY  in  the  figure  is  fixed  to  ground,  OX  being  in  the  tangent 
plane  to  the  earth  in  the  direction  of  aim.  A  second  system  denoted  by  xy  is 
also  fixed  and  in  the  XY-plane;  however,  the  x-axis  of  this  set  of  coordinates 
makes  an  angle  y  +  6  ^  with  X  .  This  total  angle  is  the  angle  of  elevation  when 
the  structure  is  loaded  with  its  dead  weight.  Finally  the  attached 

to  the  launcher  and  rotates  with  it.  This  angle  of  rotation  which  will  be  called  6 
represents  the  deviation  of  the  launcher  from  its  equilibrium  position  due  to  its 
weight. 

The  rail  RS  ,  has  a  slight  curvature.  Its  equation  in  the  system  will 
be  taken  in  the  form 


FIGURE 


5. 


where  e  Is  small  and  of  the  order  of  magnitude  of  the  angle  6  .  If  ) 

are  the  coordinates  of  the  mass  center  of  the  rocket  In  the  |ti- system,  the  angle 
of  rotation  ^  of  the  missile  relative  to  this  set  of  coordinates  Is  given  by 

If,  s  € - -  . 

To  obtain  the  equations  of  motion  for  the  system  we  first  consider  the 
free  body  diagram  of  the  missile  shown  in  Fig.  (3).  These  equations  will  be 
established  relative  to  the  -coordinates.  The  following  forces  act  on  the  rocket: 


T  (t)  -  the  thrust  force  with  misalignment  parameters  and  5  , 


W_ 


m 


€12  sin  (  12  t  +  X  )  -  shaking  force  due  to  mass  unbalance  in 

g  ■*.  2 

SOSR,  rocket, 

2 

I  12  sin  (12  t  +  X  )  -  shaking  moment  due  to  mass  unbalance  in 
xy  1 

SOSR  rocket, 

Q,  P  -  front  and  rear  shoe  loads, 

M  IQI  ,  M  I  P]  *  friction  forces  at  the  shoes. 


All  of  the  quantities  6,1  are  assumed  to  be  small  so  that  we  may 

X  A 

neglect  products  of  these  with  6  ,  ^  ,  or  e  . 

For  the  motion  of  the  mass  center  of  the  missile  we  have 
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-  /*  €  [  |P|  £•  (i^  -  a)  +  IQI  f  (e^  +  b)]  -  cos  y 
W  2 

+  R  sin  (  n  t  +  ^2  )  +  T  (  6  +  0  )  .  (2) 
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The  angular  motion  of  the  rocket  is  expressed  in  the  following  relation 


!„,(  ®  =Qb -Q€hf(e^+b)  -M  IQI  h-MiQl  €f(e^+b)b 

-Pa-P€f(e^-a)h-M  |p1  h+M|P|  €  f  -  a)  a  -  T  £2  '  T6  (a+d) 
2 

+  I  R  sin  (  R  t  +  \.  )  .  (3) 

xy  i 

These  equations  contaiii  the  unknowns  6  ,  P,  Q.  One  final  equation  is 

sufficient  for  a  unique  determination  of  these  variables.  The  moment  of  momentum 
principle  with  respect  to  the  point  O  for  the  entire  system  will  give  one  such 
equation: 
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Although  two  more  independent  equations  exist  for  the  launcher,  we  shall  not 
be  concerned  with  them  here  since  they  only  serve  to  define  the  unknown  forces 
acting  at  the  pivot  point  O. 

For  the  model  considered  we  shall  employ  Eqs.  (1),  (2),  (3)  and  (4)  to 
justify  the  approximations  which  will  be  presented  in  the  next  chapter. 


9. 


II.  APPROXIMATE  ANALYSES 

The  inaccuracies  described  in  Chapter  I  enter  into  the  analytical  descr4)tion 
of  the  system  through  non-zero  values  of  the  parameters  e,  e  ,  6,  I  ,  /i . 

A  A  Xy 

We  disregard  coupling  effects  and  examine  the  consequences  of  takiig  these 
quantities  to  be  non -zero  one  at  a  time. 

Idealized  Case  (all  parameters  zero) 

As  our  first  problem,  let  us  consider  approximating  the  response  of  the 
ideal  system:  the  one  in  which  all  of  the  inaccuracy  parameters  are  zero.  In 
this  case  Eqs.  (1)  and  (4)  become 


^  r  -  (ti  +  h)  0  1  =  T  -  sin  y  =  T  (t)  T 
g  U^m  '‘o'J  m  ^  e''o 

W 

^[0(1,+ 1  +  — )1  +^0  =  -(g  -i  ;W  cosy 

at  L  '  L  m  g  ^m  'J  ^  '^m  '’mo'  m  ' 


(5) 

(6) 


We  have  neglected  first  order  changes  in  moment  arms  for  W^  and  also  assumed 
that  p  is  very  large.  Neglecting  the  effects  of  0  in  Eq.  (5)  we  have 

^m“  ^  [(t-T)Tg(T)dT 

m  -0 

and  so  Eq.  (6)  becomes 

g^[0(l  +  ^  e^)]+„2Q,.  J^(t-T)  Tg  (T)  dT  (7) 

2 

where  I  ®  +  Ijn  «  =  p/l  .  Integrating  this  equation  with  respect  to  t , 

there  results 

.  7  7  T  g  cos  y  ^  , 

e(l  +  -T^  e^)  +  «  /0(T)dT= - 2^.^ -  /(tVTg(T)dT  (8) 

**  -'o  -'o 
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We  have  used  the  fact  that  §  (0)  »  0.  This  equation  is  essentially  an  impulse- 
momentum  relation  involving  angular  quantities.  The  first  term  is  a  change  in 
angular  momentum  reduced  by  a  factor  I,  and  the  remaining  terms  are  angular 

impulses  reduced  by  the  same  quantity.  These  impulses  represent  the  effects  of 

* 

the  spring  /3  and  the  thnist.  In  order  to  estimate  6  at  tip-off,  it  is  necessary 
to  know  the  dependence  of  6  on  time  up  to  tip-off. 

Previous  numerical  studies  along  with  field  measurements  indicate 
that  moat  realistic  systems  have  a  S-time  history  as  shown  in  Fig.  (5). 


FIGURE  5 


Since  according  to  Eq.  (7)  and  the  Initial  conditions  the  second  and  third 
derivatives  of  0  vanish  at  t  =  0,  we  shall  assume  that 


(9) 


11. 


where  6p  is  the  pitch  velocity  at  tip-off  time  tp.  Placing  Eq.  (9)  in  Eq.  (8) 
and  performing  the  obvious  reductions  we  find  that 


e  = 
P 


Tq  g  cos  y 

~n 


(t'T)  T  (T)dT 
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W  ,  (ut  > 

g  I  '  p' 


(10) 


For  systems  with  large  launch  time  (t)  may  be  replaced  by  unity  and  in  this 
case  Eq.  (10)  simplifies  to 


6 
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Tq  g  cos  y 
— FT - 


1  -I- 


(11) 


A  comparison  of  this  formula  with  known  results  indicates  that  it  is  fairly 
accurate  when  the  time  of  launch  tp  is  one-half  or  less  than  the  natural  period 
of  the  system;  that  is, 


t  -  T 
V  ~  2  n 


1 

2  w 


The  reason  for  this  can  be  obtained  by  converting  Eq.  (8), which  we  write  more 
generally  as 

9  A(t)+w^  /edT  =  F(t)  ,  (12) 

-'o 


into  an  integral  equation  for  6  and  then  solving  this  -equation  by  iteration. 
In  Eq.  (12) 


A(t)  = 


(t) 


12. 


£md  F  (t)  is  any  forcing  function.  (As  we  shall  see  later  most  of  the  inaccuracies 
present  themselves  as  different  forcing  functions. ) 

By  noting  that 

fe  (t)  dr  =  [(t-T)  6  (t)  dr  (13) 

■^0  ^  . 

if  0  (0)  =  6(0)  =  0,we  have  immediately  that  Eq.  (12)  becomes 

A  (t)  e  (t)  +w^  ^  dT  =  F  (t)  (14) 

which  is  now  an  integral  equation  in  6  .  If  now  we  change  variables  according  to 

t  t 

t  =  t  u,  —5-2  A  (t)  =  B  (u)  .  F  (t)  =  G  (u),  ^  =  a,  (1^ 

P  (2  n/  (2nr  n 

then  we  have 

-  u 

B<“)  ar 

2 

The  solution  of  this  equation  can  be  expanded  in  a  power  series  in  a  according 
to 

de  _ 

TTT  ■ 

where 

B(u) 
and 

B(u) 

In  particular  taking  (t)  =  1  we  find 
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When  a  is  sufficiently  small  the  first  tenn  in  Eq.  (18)  gives  6  (tp  )  fairly 
accurately.  Since 
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we  can  write 
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when  0/20  B  (1)  is  small.  But  this  expression  is  identical  with  Eq.  (11). 


dv 


14. 


Numerical  Solution  of  Eg.  (16) 

For  those  cases  in  which  this  two  term  approximation  is  not  accurate 
enough, additional  Iterations  become  prohibitive.  The  nature  of  Eq.  (16)  lends 
itself  readily  to  desk  calculator  computations  if  one  converts  it  into  finite 
difference  form.  To  do  this  we  first  break  up  the  interval  (0,  1)  into  n  equal 
subdivisions  l/n  units  long.  If 


d  e 

d  u 


u  =  k/n 


G.  =  G(  -  ) 
k  '  n  ' 


then  Eq.  (16)  becomes 

k-1 

0’j^  t  (  aVn^)  y  (k-m)  6*^^  =  Gj^  k  =  2 . n 

m^l 

e'o  =0  B^e’i=Gj 


(19) 


This  set  of  recursion  relations  is  readily  evaluated  on  a  desk  computer  for  6'^^ 
which  corresponds  to  the  tip-off  pitch  velocity. 

We  shall  now  show  how  to  reduce  the  effects  of  the  various  disturbing 
agencies  to  an  integral  equation  which  can  be  readi  ly  solved  either  by  iteration 
or  recourse  to  the  numerical  procedure  of  Eq.  (19). 


Thrust  Misalignment 

After  ignition  of  the  motor,  the  thrust  may  not  act  ±rough  the  mass 
center  of  the  missile.  Our  interest  is  in  evaluating  the  difference  in  pitch 
velocity  between  the  ideal  system  and  the  one  in  which  thrust  malalignment  is 
present.  If  A  ^  -  6  where  6,^  is  the  angular  response  of  the  structure 

with  thrust  misalignment  present,  then 


15. 


A  (t)j  +  <ji  A,j. 


T  6  (a  -hd  )  ^ 


(20) 


Eq.  (20)  is  the  differential  equation  satisfied  by  the  variation  In  pitch  velocity 
due  to  thrust  misalignment.  This  is  readily  converted  into  an  integral  equation 
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(21) 


When  t  /T  is  small  two  iterations  are  sufficient  for  a  solution  and  we  find  that 
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for  the  case  of  constant  thrust. 

For  those  cases  in  which  this  two  term  approximation  is  not  valid  a  solution 
of  Eq.  (21)  based  on  the  process  of  Eq.  (19)  is  in  order. 


Shaking  Forces 

When  the  missile  spins  on  the  rail  the  equation  for  the  deviation  in  pitch 
velocity  A  „  =  9  -  9  is  of  the  form 

9  S 

nr  19  I 

^[A(t)  Ag(t)J+  -  sin(n  t+\^) 

W  j 

+  — -p-  6^  n  sin  ( S2  1+^2)  (23) 

where  ^2  random  angles  between  0  and  2  ir  .  This  is  again  reduced  to 
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an  Integral  equation  in  A  „  ; 


A(t)  Ag  + 


f(t-r)  i 
JO 


Ag  dr  =  Ag  (0) 


10  .  w  2r^ 

+— - j^cos  -  cos  (  fi  t  +  Xj^)  J  +  Cj  n  jS  (t)  sin  (0  t  +X2)dT 


(24) 


The  solution  of  this  equation  by  iteration  yields  very  unwieldy  e:q)res8ions.  Further¬ 
more  because  of  the  random  character  of  Ag(0),  in  itself  is  in¬ 

determinate.  It  is  more  appropriate  to  speak  of  the  mean  value  of  A  ,  "K  .  and 
its  standard  deviation  (t  [  A  g]  .  These  can  be  obtained  by  first  noting  that 
Eq.  (24)  can  be  written  in  the  form 

A  (t)  Ag  +  j  (t-T)  Ag  dT  =  Ag  (0)  (t)  +  (t)  sin  X^  +  (t)  cos  X^ 

Jq 

+  Lg  (t)  sin  X2  +  (t)  cos  X2 

where  LQ(t)s  1.  Writing 

Ag  (i)  =  Ag  (0)  Ag^  (t)  +  Ag^  (t)  sin  X^  +^g2(*)  ^2 


+  Ag4(t)  cos  X2 


(25) 


we  find  that 


A(t)Agi 


|(t-T)  Ag^ 


dr  =  L^(t) 


(26) 


Each  of  these  can  be  solved  for  Ag^  (tp  )  by  the  process  outlined  in  Eq.  (19). 

Consequently  the  mean  value  and  standard  deviation  of  A  g  (tp  )  due  to  the 

initial  value  A  (0)  is 
8 


! 

I 

i 

I 

j 


17. 


Initial  value 


[^B<‘p>] 


initial  value 


and  for  the  parameter 


[i,(tp)]  =0, 


<'p)+^s2<'p) 


and  for  the  parameter  \2 


^2 


''[a.  Up)] 
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83  '  p  '  s4  '  p  ' 


Consequently  the  expected  mean  value  of  ^  ^  (tp  )  is 

IJiT)  -  pT^]  A^„(.p) 


and  the  standard  deviation  is 


"  [  ^p(‘p  »]  -'AA'  ^L<V^^[‘^b1<V^^b2<V^^b3<V^^s4<‘p>]- 

In  computing  these  we  have  assumed  that  \2*  independent. 


F  riction 


In  treating  friction  we  assume  that  is  given  by  the  expression 


m  -Jo 
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and  write 
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where  N  is  taken  as 


So  long  as  N  a*  0  we  may  write 


e. 


A(t)+n  -f- 


’m 


2  ' 
1“  (t-T) +i^— r—  — 
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t  -  f  +h)  COS-y  t  . 


An  approximate  expression  for  8^  follows  by  iteration.  When  this  is  done 
there  results 
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A  more  accurate  solution  can  be  obtained  by  writing  Eq.  (28)  as 
■t 

c(t)6j+  /  D(t.T)e,(T)dT=  E(t) 


and  solving  Eq.  (29)  numerically.  By  converting  to  the  variable  t=  tp  u 
have 

r 

L  (u)  e’j  (u)  +  /  M  (u,  ' )  ^  d  V  =  N  (u) 

■Jo 

Splitting  the  range  (0, 1)  into  n  equal  subintervals  there  results 


0^(7)  dr 
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Rail  Curvature 


As  a  final  item  we  consider  rail  curvature.  In  this  case  all  we  present 
is  an  integral  equation.  Any  solution  of  this  equation  hinges  on  an  assumed  form 
for  rail  curvature.  If  Aj^  denotes  the  difference  ®  where  0^^  is  the 
angular  response  in  the  presence  of  rail  curvature  Eq.  (4)  reduces  to 
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(30) 


It  is  worth  noticing  that  rail  curvature  introduces  a  thrust  misalignment  in  a 

sense.  The  appropriate  Integral  equation  is 
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This  is.  again  of  the  form  which  Bq.  ( 19)  solves  numerically. 
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